Measurement of the Muon Decay Parameter 8 
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The muon decay parameter 5 has been measured by the TWIST collaboration. We find 5 = 
0.74964 ± 0.00066(stat.) ± 0.00112(syst.), consistent with the Standard Model value of 3/4. This 
result implies that the product P^^ of the muon polarization in pion decay, P M , and the muon 
decay parameter £ falls within the 90% confidence interval 0.9960 < P M £ < £ < 1.0040. It also has 
implications for left-right-symmetric and other extensions of the Standard Model. 

PACS numbers: 13.35.Bv, 14.60.Ef, 12.60.Cn 



The TWIST spectrometer was designed to mea- 
sure a broad range of the normal muon decay spectrum, 
jjL + — > e + v e V^, allowing the simultaneous extraction of 
the spectrum shape parameters. Assuming the weak in- 
teraction is local and invariant under the Lorentz group, 
the effective four fcrmion muon decay matrix element can 
be written in terms of helicity-prescrving amplitudes: 
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where the g]^ specify the scalar, vector, and tensor cou- 
plings between ^i-handed muons and e-handed electrons 
|2(. In this form, the Standard Model implies g\ L = 1 
and all other coupling constants are zero. 

The differential decay spectrum Q of the e + emitted 
in the decay of polarized p + is provided in terms of four 
parameters, p, 5, r\, and £, commonly referred to as the 
Michel parameters, which are bilinear combinations of 
the coupling constants. In the limit where the electron 
and neutrino masses as well as radiative corrections are 
neglected, this spectrum is given by: 

^M^o) * 3 ( 1 -*) + ^ 4 *- 3 > 

+ P^cos6[l-x+^5{4x-3)}, (2) 

where 9 is the angle between the muon polarization and 
the outgoing electron direction, x — E e /E max , and is 



the muon polarization. The fourth parameter, 77, appears 
in the isotropic term when the electron mass is included 
in the analysis. In the Standard Model, the Michel pa- 
rameters take on precise values. 

The parameter £ expresses the level of parity violation 
in muon decay, while <5 parametrizes its momentum de- 
pendence. Recently, TWIST reported a new measure- 
ment of p In this paper we report a new mea- 
surement of 6. The currently accepted value of 6 = 
0.7486 ± 0.0026 ± 0.0028 agrees with the Standard 
Model expectation of 3/4. Some Standard Model exten- 
sions require deviations from pure V — A coupling that 
can alter 6. Some of these models involve right-handed 
interactions. The positive definite quantity, 
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can serve to set a model independent limit on any muon 
right-handed couplings 0, 01 . A recent review of muon 
decay is presented in . 

Highly polarized surface muons are delivered to the 
TWIST spectrometer Q from the M13 channel at TRI- 
UMF. The spectrometer consists of a detector made up of 
56 very thin high precision chamber planes, all mounted 
perpendicularly to a solenoidal 2 T magnetic field. The 
muons enter this array of chambers through a 195 pm 
scintillator that acts as the event trigger. More than 80% 
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of the mucins come to rest in the central stopping target, 
which also acts as the cathode plane for the Multi Wire 
Proportional Chambers (MWPC) on either side. The de- 
cay positrons spiral through the chambers producing hits 
on the wires that are recorded by time to digital convert- 
ers. These helical tracks are later analyzed to determine 
precisely the positron energy and angle. The observed 
momentum resolution is 100 keV/c |4|. The cos 9 reso- 
lution derived from Monte Carlo (MC) is about 0.005. 
The reconstruction is similar to except for some de- 
tails discussed below. 

TWIST determines the Michel parameters by fitting 
two-dimensional histograms of reconstructed experimen- 
tal decay positron momenta and angles with histograms 
of reconstructed Monte Carlo data. This approach has 
several advantages. First, spectrum distortions intro- 
duced by the event reconstruction largely cancel because 
MC and experimental data are analyzed identically. Sec- 
ond, because the MC simulates the detector response 
well, no explicit corrections of the result are required. 
Third, a blind analysis of the result is straightforward. 
It is implemented by utilizing hidden Michel parameters 
Ph,$h, and £h to generate the theoretical decays. The 
decay rate can be written as 
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since the decay spectrum is linear in the shape param- 
eters. The sum of MC spectra is fit to the data spec- 
trum by adjusting the AA. 5 is extracted as (^hSh + 
A (££))/(£» + A£). Since the hidden parameters were al- 
lowed to only deviate from their Standard model values 
by no more than 0.03 it was sufficient for the extrac- 
tion of systematic uncertainties to assume that they had 
their Standard model values during the blind stage of the 
analysis. The MC spectra were generated including full 
0(a) radiative corrections with exact electron mass de- 
pendence, leading and next-to-leading logarithmic terms 
of 0(a 2 ), leading logarithmic terms of 0(a 3 ), corrections 
for soft pairs, virtual pairs, and an ad-hoc exponentiation 
0. Because TWIST at the present stage could not pro- 
vide an improved measurement of eta, we set it, for MC 
spectra production, to its current highest precision value 
of —0.007 in order to constrain other parameters bet- 
ter. The uncertainty of 0.013 on the accepted value of rj 
gives a negligible uncertainty on the final value of 5. 

The TWIST simulation model is based on GEANT 
3.21 ^} with the chamber response based on 
GARFIELD It contains virtually all the compo- 

nents of the spectrometer with which a muon or a decay 
positron could interact. The output exactly mimics the 
binary files generated by the data acquisition system. 

Factors that influence the momentum and angle deter- 
mination must be well simulated in the MC, so special 
runs were taken specifically to address the accuracy of the 



| 400 
O 300 



200 



100 




-2 -1.5 -1 -0.5 0.5 1 13 2 
Change in Momentum (MeV/c) 



CO 



o 
O 



200 


-AG 


1 

tlLt 1 • Data 


150 






100 










V 


50 


- ^ 








inHfilii 



-0.04 -0.02 0.02 0.04 
Change in Angle (radians) 



FIG. 1: (color online) The difference between downstream 
and upstream tracks, for both data and MC, resulting in: 
Top, the positron momentum change in the central stopping 
target, Bottom, A9 for a positron that passed through the 
central stopping target. The MC results were normalized to 
those of the data for the purposes of this figure. 



simulation of energy loss and multiple scattering. Muons 
were stopped in the extreme upstream wire chambers in 
both the experiment and in the MC simulation. The 
decay positrons were tracked through the upstream half 
and separately through the downstream half of the spec- 
trometer. Differences in momentum and angle were his- 
togrammed on a track by track basis. Figure ^ presents, 
for both data and MC, the changes in momentum and an- 
gle that occur primarily at the central stopping target. 
The widths of the peaks in this figure do not represent the 
experimental resolution for a number of reasons. First, 
because the same track is being reconstructed twice with 
finite resolution the differences in the measured values 
can be either positive or negative. Second, the particle 
sees approximately twice the thickness of materials. As 
well, the tracking in the upstream region is in the op- 
posite direction for which the code is optimized and in 
this region the track does not see as many planes due 
to the distribution of the muon stops in the upstream 
planes. The MC nonetheless reproduces the data very 
well. The Ap distribution mean(RMS) for the data and 
MC arc -0.17(0.41) McV/c and -0.17(0.39) MeV/c, re- 
spectively. The A9 mcan(RMS) for the data and MC 
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are -0.95(17.0) and -0.37(18.0) milliradians respectively. 
The small differences are within the uncertainties associ- 
ated with positron interactions and target thickness. 

The result for 5 presented here employed a sample con- 
sisting of 6 x 10 9 events recorded in Fall, 2002. This data 
sample is comprised of the same 16 data sets used for our 
extraction of p 0]. Many of these data sets were taken 
under conditions chosen to establish the sensitivity of the 
detector to systematic effects. Four of the data sets, sets 
A and B taken at 2.00 T six weeks apart and two other 
sets, one taken at 1.96 T and one at 2.04 T, were ana- 
lyzed and fit to their corresponding MC samples to derive 
the value of 8. Our p determination also utilized a cloud 
muon sample 0. The low polarization of that data set 
leads to low sensitivity and the potential for substantially 
increased systematic effects in the extraction of 8. Rather 
than perform a complete additional systematics study for 
a data set that would contribute little weight, we chose 
to use the cloud muon sample only as a consistency check 
on our final result for 8. 

There are several differences between our previous 
analyses for p Q and for the 8 result presented here. 
After the first analysis was completed and the hidden pa- 
rameters were disclosed, an a priori defined consistency 
check was carried out. Muon decay parameters deter- 
mined from set B were used to generate a new Monte 
Carlo spectrum, which was used to perform another fit 
to that set. This fit was expected to yield deviations of 
all parameters consistent with zero, but it failed for 8. 
It was determined then that there was a flaw in the way 
the polarization-dependent radiative corrections were im- 
plemented in the Monte Carlo event generator. Since p 
is essentially decoupled from the asymmetry parameters 
this flaw had no impact on the value of p. After the first 
analysis, we were no longer "blind" to the value of p. 
However the flaw did introduce a systematic uncertainty 
in the value of 8 inferred from the data. This systematic 
uncertainty depended on the difference in absolute po- 
larization (not P M £) between the actual muon decay data 
and the Monte Carlo decay events. Estimates indicated 
that the systematic effect was < 0.001. Unfortunately, 
a precise value was impossible to determine a priori, as 
only the product P M £ is measurable from the data, rather 
than P M alone. 

A new analysis with a corrected event generator was 
therefore undertaken for 8 through the generation, with 
a new set of hidden Michel parameters, of a completely 
new set of Monte Carlo events. The four data sets de- 
scribed above were also reanalyzed with improved align- 
ment calibrations. Finally, the track-selection algorithm 
was improved by merging those used for the p analysis. 
Thus the fits employed for the 8 extraction are completely 
distinct from those in p|. 

Figure [21 shows the decay positron angular distribu- 
tions for representative momentum bins. Equation @ 
indicates the angular distributions follow a 1 + A(p) cos 8 
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FIG. 2: (color online) Decay positron angular distributions 
from set B (solid curves) and the corresponding best fit distri- 
butions within the fiducial region (dashed curves) for selected 
momentum bins. Z = ir — 0. 



shape, where by convention the asymmetry, A(p), is pos- 
itive when positrons are emitted preferentially along the 
muon polarization axis. Figure |3f a) shows the observed 
muon decay asymmetry as a function of momentum for 
set B. The asymmetry provides a compact representation 
of the angular distributions. However, extracting 8 from 
A(p) involves a significant correlation between 8 and p 
In contrast, we extract 8 from a simultaneous fit of 
the full experimental momentum-angle distribution illus- 
trated in Fig. as described above, which leads to a 
negligible correlation between 8 and p. Fits to upstream 
minus downstream distributions, which are essentially in- 
dependent of p and 77, gave nearly identical results for 8. 

The fiducial region adopted for this analysis requires 
p < 50 MeV/c, \p z \ > 13.7 MeV/c, p T < 38.5 MeV/c, 
and 0.50 < | cos6*| < 0.84. The fiducial cuts, while inten- 
tionally chosen to be conservative, are related to physical 
limitations of the TWIST detector. The 50 McV momen- 
tum cut rejects events that are near the region utilized 
in the end point fits Q. It is also important to avoid 
the region very close to the end point to minimize the 
sensitivity of the Michel parameter fits to the momen- 
tum resolution. The longitudinal momentum constraint 
eliminates events with wavelengths that match a 12.4 cm 
periodicity in the wire chamber construction. The trans- 
verse momentum constraint insures that all decays are 
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well confined within the wire chamber volume. The an- 
gular constraint removes events at large cos 9 that have 
worse resolution and events at small cos 9 that experience 
larger energy loss and multiple scattering. These limits 
were frozen early in the analysis. Prior to opening the 
"black box" a study of how the results changed as each of 
the fiducial boundaries was moved found the sensitivities 
to be very weak. 

Figure|21shows the results of the best fit to set B within 
the fiducial region for the selected momentum bins. Set 
B is one of the statistically larger sets and is typical of 
all the sets. Figure shows the measured muon decay 
asymmetry in panels (a) and (b) while panel (c) presents 
the difference between the measured asymmetry and the 
asymmetry calculated from the best fit MC spectrum for 
events within the fiducial region. Panel (d) shows the 
difference between the asymmetry within the fiducial as 
reconstructed and as thrown for the MC, illustrating that 
the distortion of the asymmetry by the TWIST detector 
is small and essentially momentum independent. 

The graphite coated Mylar stopping target resulted in 
a time dependence of the muon polarization. P t ,. which 
prevented the simultaneous determination of a value for 
P M £ from this data sample. (P M ) ~ —0.89 at the time of 
decay for the data sets analyzed here. The graphite coat- 
ing on the Mylar target was necessary since the target, 
also serving as a cathode foil for the two central MWPC 
chambers, required a conductive surface. Details regard- 
ing the target can be found in Knowing the precise 
polarization is not important for extraction of 6, thus 
the Mylar target was considered adequate for the cur- 
rent measurement despite the possibility of depolarizing 
interactions. 

Systematics were studied by employing the fitting tech- 
nique described above to fit experimental data samples 
taken with a systematic parameter set at an exaggerated 
level to data taken under ideal conditions. This expresses 
the changes in the spectrum shape caused by the system- 
atic effect in terms of changes in the Michel parameters. 
Other systematic sensitivities were determined by ana- 
lyzing a data or MC sample with a systematic parameter 
offset from its nominal value and fitting to the same sam- 
ple analyzed with this parameter at its nominal value. 
For the current analysis the largest uncertainties are for 
the detector alignment, for the simulation of positron in- 
teractions, and for the chamber response, in particular 
the time dependent effects due to gas density changes 
and to the variability of the cathode foil positions 
The latter parameters were monitored throughout the 
data accumulation periods and average values were used 
in the analysis. Uncertainties due to the detector align- 
ment were established by analysis of data and genera- 
tion of MC with purposely misaligned chambers. Up- 
per limits for the positron interaction uncertainties were 
derived from studies of the data for muons stopped far 
upstream and from MC histograms that demonstrated 
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FIG. 3: (color online) (a) The observed muon decay asym- 
metry from set B for all events within 0.50 < | cos#| < 0.84. 
(b) The same quantity for those events that fall within the 
fiducial region, (c) The difference between the data in panel 
(b) and the best fit MC spectrum, (d) The difference between 
the asymmetry calculated from the reconstructed MC events 
and as thrown by the MC, which illustrates the spectrometer 
response. 



the distortion of the momentum spectrum due to hard 
interactions. Other important systematic uncertainties 
for S are the stopping target thickness and the momen- 
tum calibration. The target thickness issue was studied 
by varying the thickness of the graphite coating in MC. 
The results of these studies for the parameter 8 are pre- 
sented in Tables H] and [D] The value of S for the cloud 
muon sample is 0.75245 ± 0.00526(stat.), consistent with 
the results in Table I. The average value of p from the 
present fits is 0.75044, consistent with the blind analysis 
result in Rcf. 

The effects of chamber response, momentum calibra- 
tion and muon beam stability, which have time depen- 
dent components, are treated as data set-dependent ef- 
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TABLE I: Results for 5. Each fit has 1887 degrees of free- 
dom. Statistical and set-dependent systematic uncertainties 
are shown. 



Data Set 


,5 


x 2 


Set A 


0.75087 ± 0.00156 ± 0.00073 


1924 


Set B 


0.74979 ± 0.00124 ± 0.00055 


1880 


1.96 T 


0.74918 ± 0.00124 ± 0.00069 


1987 


2.04 T 


0.74908 ± 0.00132 ± 0.00065 


1947 



TABLE II: Contributions to the systematic uncertainty for 
6. Average values are denoted by (ave), which are considered 
set-dependent when performing the weighted average of data 
sets. 



Effect 



Uncertainty 



Spectrometer alignment 
Chamber response(ave) 
Positron interactions 
Stopping target thickness 
Momentum calibration(ave) 
Muon beam stability(ave) 
Theoretical radiative corrections[9j 
Upstream/Downstream efficiencies 



±0.00061 
±0.00056 
±0.00055 
±0.00037 
±0.00029 
±0.00010 
±0.00010 
±0.00004 



fects with the avcrage(ave) over the four sets used in the 
S evaluation appearing in Table ITT1 

We find 5 = 0.74964 ± 0.00066(stat.) ± 0.00T12(syst.), 
consistent with the Standard Model expectation of 3/4. 
The central value for S was calculated as a weighted av- 
erage using a quadratic sum of the statistical and set- 
dependent uncertainties for the weights. The final sys- 
tematic uncertainty is a quadratic sum of set independent 
and average values of the set-dependent systematics. Us- 
ing this result, our new value for p the previous mea- 
surement of P^S/p ^2], an d the constraint Q R > 0, it is 
possible to establish new 90% confidence interval limits, 
0.9960 < P^ < £ < 1.0040, consistent with the Standard 
Model value of 1. This result is more restrictive than the 
current best measurements for muons from pion and kaon 
decays 0, ^| . In addition, from these same results one 
finds that Q R < 0.00184 with 90% confidence. This may 
be combined with Eq. (PJ to find new 90% confidence 
limits on interactions that couple right-handed muons 
to left-handed electrons: \g[ R \ < 0.086, \gl R \ < 0.043, 
and \g1 R \ < 0.025. The lower limit, 0.9960 < P^ can 
be used to determine a new limit on the mass of the 
possible right-handed boson, Wr, improving the exist- 
ing lower limit of 406 GeV/c 2 (402 GeV/c 2 with modern 
M Wl = 80.423 GeV/c 2 ) from to 420 GeV/c 2 under 
the assumption of pseudo manifest left-right symmetry. 
For nonmanifest left-right symmetric models the limit 
is Mw r 9l/ gs, > 380GeV/c 2 , where <jx and gn are the 



coupling constants 0|. The value of 5 is sensitive to a 
proposed nonlocal interaction [lfi | that would be repre- 
sented by a new parameter k. A limit for k may be es- 
timated from our 90% confidence lower limit for 5 using 
the relation 8 = 3/4(1 - 6k 2 ). This results in k < 0.024, 
which compares with k = 0.013 [lfi | hinted at by it decay 
experiments. 
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